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Abstract. In this note we study, for a random lattice L of large dimension n, the 
supremum of the real parts of the zeros of the Epstein zeta function E n (L, s) and 
prove that this random variable has a limit distribution, which we give explicitly. 
This limit distribution is studied in some detail; in particular we give an explicit 
formula for its distribution function. 



1. Introduction 

Let X n denote the space of all n-dimensional lattices L C M n of covolume one. 
For L € X n and Jfts > the Epstein zeta function is defined by 



1.1) E n (L,s) = ^\v\ 

veL 



where ' denotes that the zero vector should be omitted. E n (L,s) has an analytic 
continuation to C except for a simple pole at s = f| with residue tt^T(^) . Fur- 
thermore, E n (L,s) satisfies the functional equation 

(1.2) F n (L,s) = F n {L*,%-s), 

where F n (L, s) := ir~ s T(s)E n (L, s) and L* is the dual lattice of L. 

The Epstein zeta function is in many ways analogous to the Riemann zeta function. 
In particular we have the relation 

E 1 (Z,s) = 2((2s). 

Because of this analogy and for other related reasons, many studies have been made 
regarding the location of the zeros of E n (L, s). From (|1.2|) it is clear that E n (L, s) has 
a "trivial" zero at each point s = —1, —2, —3, . . ., just like C(2s), and the remaining 
nontrivial zeros of E n (L, s) are in bijective correspondence with the nontrivial zeros 
of E n (L* , s) under the map s i— > § — s. However, the Riemann hypothesis for E n (L, s) 
generally fails: E n (L,s) typically has many nontrivial zeros which do not lie on the 
critical line = f . Cf. 0, [I], [23], [29], [30], [H], M- 

We denote by Nl(T) the number of nontrivial zeros (counting multiplicity) of 
E n (L,s) with \$ss\ < T. Then Nl(T) satisfies the following Riemann-von Mangoldt 
type asymptotics 



2T T 

1.3) N L (T) = — log — +O h {\ogT) asT^oc, 

7r TTem{L)m{L* ) 
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where m(L) is the length of the shortest non-zero vector in L. 

From the point of view of number theory, the most interesting choices of L are 
those for which the Gram matrix for some (and thus any) Z-basis of L is proportional 
to an integer matrix. We call these lattices rational. In particular when n = 2 many 
results have been obtained regarding the zeros of E%{L, s) for rational L correspond- 
ing to integer quadratic forms with a fundamental discriminant. It was conjectured 
by H. L. Montgomery that in this case asymptotically 100% of the nontrivial zeros 
of E2(L, s) lie along the critical line Jis = \. This was proved conditionally, assum- 
ing the Generalized Riemann Hypothesis and a weak form of well-spacing for the 
zeros of L-functions attached to ideal class characters, by Bombieri and Hejhal in 
[6]. Furthermore, Selberg has proved unconditionally, in still unpublished work (cf. 
[TTJ p. 553] and [5l pp. 225-227]) that a positive proportion of the zeros do lie on the 
critical line. For related results, see also [2] and [TTJ. 

Our main object of study in the present paper is the supremum of the real parts 
of the zeros of E n (L, s), i.e. 

a L : = sup{3fy? : E n (L,p) = 0}. 

In other words, o~l gives the precise location of the zero- free right half-plane of 
E n (L,s). One easily shows that o~l exists and is finite for any given L £ X n ; 
furthermore &l > § always holds (cf., e.g., [24, p. 693 and Thm. 1]). Of course, 
= ol* = j is equivalent with the Riemann hypothesis for E n (L,s). Note that 
o~l is lower semicontinuous (and hence Borel measurable), since any zero s = sq of 
E n (Lo,s) gives rise to a nearby zero for all E n (L,s) with L in a sufficiently small 
neighborhood of Lq (as follows from a standard application of Rouche's theorem 
using the formula [2TJ (23)] for ir~ s T(s)E n (L, s)). We also remark that o~l takes 
arbitrarily large values for any given n > 2 (cf . Remark in Section [2D . 

For any Dirichlet series f(s) = YlJLi e -? s with exponents Ai < A2 < . . . whose 
pairwise differences do not satisfy any non-trivial linear relation over Q, the supre- 
mum of the real parts of the zeros of f(s) equals the unique number a for which 
e~ Al<J = X^2 e ~ AjCT 5 c f- Lemma Q] below. This independence condition never holds 
for E n (L, s) (e.g. since L contains both 2v and 4i; for any v G L). However, we have 

(1.4) E n (L, s) = 2((2s) \v\~ 2s (3?« > f ), 

where L denotes a set containing one representative from each pair {v, —v} of prim- 
itive vectors in L; and it turns out that the Dirichlet series X^gL IH~ 2s satisfies 
the independence condition for fi n -almost every lattice L G X n , where \i n is Siegel's 
measure ([22]) on X n ; see Lemma [2l From this we conclude (cf. Section [2|): 

Proposition 1. Let n > 2. For almost every L € X n , o~l equals the unique number 
a > which satisfies 2m(L)~ 2a = ^((2a)^ 1 E n (L,a). It follows that for almost 
every L £ X n , E n (L,s) has infinitely many zeros with > \. 

In particular, for small n the formula in Proposition[T]makes it possible to compute 
a 1, numerically for a given generic L £ X n . We stress, however, that for a lattice L 
such that Y1 V i=l \ v \~ 2s does not satisfy the linear independence condition (e.g. any 
rational L, cf. Remark U] in Section [2]), the computation of aj_, is in general not an 
easy task. We mention that Bombieri and Mueller in [7] have shown how to calculate 
a 1, explicitly for certain examples of rational lattices L £ X<i (with o~l > 1), where 



ON THE ZERO-FREE HALF-PLANE OF A RANDOM EPSTEIN ZETA FUNCTION 



3 



they also obtained bounds on the asymptotic rate of approach of the zeros of Ei{L, s) 
to the line ?fts = gl- See also [5] for a related investigation of the supremum of the 
real parts of the zeros of certain other Dirichlet series. 

Our main result concerns the distribution of gl for a random lattice L in large 
dimension n. The random element L E X n will always be chosen according to Siegel's 
measure fi n , normalized to be a probability measure. The present study is motivated 
by recent investigations [27| of the value distribution of E n (L,s) for > § and a 
/x n -random lattice L of large dimension n, where the following result is established: 
Let V n denote the volume of the n-dimensional unit ball. Let V be a Poisson process 
on the positive real line with intensity \ and let T\ , T2 , T3 , . . . denote the points of 
V ordered so that < T x < T 2 < T 3 < ■ ■ ■ . Then, for any fixed c G C with Kc> ±, 

00 

(1.5) V~ 2c E n (; cn) A T(c) := 2 Tf 2c as n -)■ 00, 

3=1 

i.e. the random variable V~ 2c E n (-,cn) converges in distribution to T(c). 

The proof of (|1.5p is built on a result [26] which provides the connection between 
the lengths of lattice vectors appearing in the formula (jl.ip and the points of the 
Poisson process V . Since this result is an important ingredient also in the present 
investigation we recall it here. Given a lattice L € X n , we order its non-zero vectors 
by increasing lengths as ±i>i, ±1*2, i^3i • • -i se t £j = \vj\ (thus < l\ < 1% < . . .), 
and define 

V,-(L) :=F n ^, 

so that Vj (L) is the volume of an n-dimensional ball of radius lj . The main result in 
PU] states that, as n — > oo, the volumes {Vj{L)}JL l determined by a random lattice 
L G X n converges in distribution to the points {Tj}°5 =1 of the Poisson process V on 
the positive real line with constant intensity \. 

In view of the last two paragraphs, together with Proposition [T] and the fact that 
C(2o") — > 1 as g — > cxd, it seems reasonable to expect that as n — > 00, u~ 1 gl should 
tend in distribution to 

00 

(1.6) a {Tj} := the unique g > \ satisfying T-f 2<T = Tr 2 °" 

(We will show in Section [3] that C{Tj} is a well-defined random variable.) Our first 
theorem states that this is indeed the case. 

Theorem 1. If L is taken at random in X n according to /%, then 

-1 d , 
n gl — > a {Tj} as n ~ > °°- 

Our second theorem gives an explicit formula for the distribution function of 0"{t }• 
Recall that the lower incomplete gamma function 7(5, z) is defined by 

(1.7) 7(s,z) := f u s ~ l e- u du 

Jo 

for s,z£C with IRs > 0. In order to make 7(5, z) single-valued, we will always keep 
z£C \ M<o (in fact, we will only need to use z with 3lz > 0), and choose a path of 
integration in (jl.7p which stays inside this cut plane. We agree that | argn| < it for 



4 



ANDREAS STROMBERGSSON AND ANDERS SODERGREN 




Figure 1. A graph of the density function of C{t }• It was computed 
using the method described in Appendix [X] (see also |25} numden- 
sity.mpl]). 

all u £ C\ M<o- Now the function 7(5, z) is extended to all s £ C \ Z<o, z £ C\ M<o 
through the recursion formula 

7(s + 1, z) + z s e~ z 



1.8) 



7(s,z) 



Theorem 2. For any c > i, we /ia?;e 



Probia 



<c) 



1 2c 

2 + V 



'0 V 7 ^'"^/ 

TTie integral in the right hand side is absolutely convergent. 



dy. 



Corollary 1. The random variable 0"{t } has a continuous density function /(c) 
given explicitly in (|5.ip below. It satisfies 



2-K 1 (c-±y + 0((c 



as c 



(1.9) 
and 

(1.10) /(c) = K 2 c~' i + O I r j «.s c — x . 

where K\ = 39.47841 . . . and K2 = 0.822467 . . . are positive real numbers given 
explicitly below in (|5.20|) and (I5.14p . respectively (cf. also (|5.3p . (|5.13|) . (|5.7p and 

to;- 

In Appendix |Aj we also give formulas for the distribution and density functions of 
C{T } obtained through the residue theorem, and discuss numerical evaluation. See 
Figure 1 for a graph of the probability density function of o~{t } generated using the 
formulas in Appendix [A] (cf. \25\ numdensity.mpl]). 

We conclude by remarking that, as is rather clear from the previous discussion, 
the random variable o"{t,,} can a l so t> e interpreted as the supremum of the real parts 
of the zeros of the random Dirichlet series f{Tj}i s ) = YlJLi^f 23 ■ Indeed, by the 
strong law of large numbers the series /{t 3 }( s ) nas > with probability one, abscissa 
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of absolute convergence ao = | and satisfies lim cr _ > i + f{T}(°~) > 2T X ; also with 
probability one the numbers 2(logTj — logTi), j = 2, 3, . . ., are linearly independent 
over Q; this means that Lemma [T] below applies almost surely and the claim follows. 
Hence Theorem [2] and Corollary Q] describe explicitly the distribution of the location 
of the zero-free right half- plane of /{t}(s). 

There exists a vast literature on random Dirichlet series; however, we are not 
aware of many results pertaining to their zeros. Cf., however, Edelman and Kostlan 
[9l §§3.2.5, 8.2], regarding the zeros of the random Dirichlet series Ylri=i a nn~ s , 
where a n are independent standard normal random variables. 

1.1. Acknowledgments. We are grateful to Daniel Fiorilli and Svante Janson for 
inspiring discussions and helpful remarks. The second author thanks the Institute 
for Advanced Study for providing excellent working conditions. 

2. Proof of Proposition [T] 

Lemma 1. Consider any Dirichlet series f(s) = Yl'jLi with real exponents 

Ax < A2 < • • • and abscissa of absolute convergence ao < 00. Assume lim o ._ j>o .+ f(o~) > 
2 e -Aio-o_ Then the equation f(o~) = 2e~ Xlcr has exactly one real root a = cry > ao- If 
furthermore all the differences Xj — Ai for j = 2, 3, . . . are linearly independent over 
Q, then o~f equals the supremum of the real parts of the zeros of f{s), and the function 
f(s) has infinitely many zeros in any strip o~\ < < 02 with ao < o\ < 02 < 07. 

Remark 1 . The linear independence condition of the lemma is equivalent to the state- 
ment that if c\ , C2 , . . . are any integers all but finitely many vanishing and satisfying 
Ylri=i c n = and Ylri=i c "An = 0, then c\ = c<i = . . . = 0. This is also equivalent 
to the statement that the pairwise differences among Ax, A2, • • • do not satisfy any 
non-trivial linear relation over Q, i.e. if Cjk for 1 < j < k are integers all but finitely 

many vanishing and satisfying Y,j<k c jk(^j ~ h) = 0, then Efelj+i c jk = E£=x c kj 
for all j > 1. 

Proof. The fact that the equation f(o~) = 2e _AlCT has exactly one real root a = aj > 
o"o follows since the function a h4 e Xl<T f(a) for a > is strictly decreasing, tends to 
1 as o~ — y 00, and by assumption tends to a limit which is greater than 2 as a — > o~q. 
For any s with 3fo > a f we have | £~ 2 e ~ Aj 1 < Y J f=2 e ~ X ^ S < e~ AlRs = |e~ AlS |, 
and thus /(s) 7^ 0. Hence it now only remains to prove that /(s) has infinitely many 
zeros in any strip o~\ < < o"2 with ao < o~\ < 02 < o~f. Assume the contrary; then 
there even exist some 01, 02 with ao < 01 < 02 < 0/ such that f(s) has no zero in 
the strip o\ < 5?s < 02- By basic facts in complex analysis, this implies that 

(2.1) inf|/(<r + it)|>0 

£<EK 

for any fixed a G (01,02) (cf. H §4 (Hilfssatz 3)] and [El §3 (Hilfssatz 3)]). 

On the other hand, for any a £ (01,02) we have YIJL2 e^ x ~^ a > 1, and hence 

there exist C2, . . . 6 C satisfying |C 2 | = |Cs| = • • • = 1 and ^^2 e^ 1_A ^ CT Cj = -1- 
It follows from Kronecker's density theorem (cf., e.g., [151 Prop. 1.5.1]), using our 
linear independence assumption, that for any given J G Z>2 and e > there exists 
tGl such that | e ( A i- A j)** _^-| < e for all j G {2, . . . , J}. Applying this with J -> oo 
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and e — > 0, we conclude that 

oo 

inf \f(<r + it)\ = e~ Xl(7 mi 1 + V e (*i-^) =0. 

i=2 

This contradicts ()2.1|) . and hence the lemma is proved. □ 
Recall the definition of the set L from just below equation (|1.4p . We now prove: 

Lemma 2. For each n > 2 and almost all L G X n the following holds: The vector 
lengths \v\ for v G L are all distinct, and the pairwise differences of their logarithms 
do not satisfy any non-trivial linear relation over Q. 

Proof We realize X n as the homogeneous space SL(n, Z)\SL(n, R), where SL(n, Z)g 
corresponds to the lattice Z n g C R n . Note that /i n is the unique probability measure 
on X n induced from a Haar measure on SL(n,R). We will let fj, n denote also the 
corresponding Haar measure on SL(n,R). Now the statement of the lemma is equiv- 
alent to the following: For almost every matrix M G SL(n,R), any finite sequence 
of primitive vectors U\, ... , ttjy G Z n (N > 2) with Uj ^ iii^ for j ^ k, and any 
bi, . . . , i>AT G Z \ {0} with X^Li ^' = 0' we have 

iV 

(2.2) J^ftjlog l^-Ml / 0. 

i=i 

Since there are only countably many possible 2A^-tuples . . . , itjv, b±, . . . , fr/v) 
it suffices to prove that for each fixed choice of (ttj, . . . , mat, 6i, . . . , &at) the set of 
M G SL(n,R) satisfying (g2D has full measure in SL(n,R). We note that is 

equivalent to rij=i 



M\ 2b > 


+ 1, i.e. 




N 




AT 


n 


1 )i/rl 26 i 


n 


i=i 






(6j>0) 







(2.3) II \ Uj Mr j - II |wjM| 2|6jl ^0. 



Hence, by the explicit formula for the measure [i n on SL(n, R) in terms of the matrix 
entries (cf., e.g., [32]) and the fact that the left-hand side of (|2.3j) is homogeneous in 
M (since X^j^i bj = 0), we find that it is enough to prove that for any fixed choice 
of (ui, . . . , ujy, bi, . . . , 6jv) as above, the relation (|2.3|) holds for Lebesgue almost all 
matrices M G Mat n , n (M) = R n2 . 

Note that each factor |ujM| 2 in ()2.3p is a real polynomial in the n 2 matrix entries 
of M. Our conditions on u\, . . . , imply in particular that U\ is not proportional 
to any of U2, ■ ■ ■ , wjv, and thus the set S of vectors in R™ which are orthogonal to 
ui but not orthogonal to any of v>2, ■ ■ ■ ,wjv is non-empty. Now, if we take any 
M G Mat nin (R) all of whose column vectors lie in S, we note that the left-hand 
side of (|2.3|) is non-zero. This proves that the left-hand side of (|2.3|) is a real, non- 
zero polynomial in the n 2 matrix entries of M. Hence the condition (|2.3p is indeed 
fulfilled for Lebesgue almost all M G Mat n>n (R), and the lemma is proved. □ 

Proof of Proposition^ Take any L G X n such that the vector lengths \v\ for v G L 
are all distinct, and the pairwise differences of their logarithms do not satisfy any 
non-trivial linear relation over Q. By Lemma[2]this holds for almost every L. Having 
fixed any such L, we consider the Dirichlet series f(s) = ^„ e £ M -2 *- The abscissa 
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of (absolute) convergence for f(s) is uq = ^, and we have lim cr _ > n+ f(a) = oo; cf. 

()1.4p . Hence, by Lemma [JJ the supremum of the real parts of the zeros of f(s) 
equals the unique real root aj > ^ of the equation f(a) = 2|i>i|~ 2<T , where v\ is the 

shortest vector in L; in fact /(s) has infinitely many zeros in any strip a\ < < a 2 
with I < o\ < (72 < (Tf. Using \v\\ = m{L) and f j 1 . 4 [) . we see that the equation 
for Of may equivalently be expressed as 2m{L)~ 2a = ^((2a)~ 1 E n (L,a). Finally, 
using (|l,4p and the fact that £(s) does not have any zeros when $ts > 1, we see 
that E n (L,s) has exactly the same zeros (also counting multiplicity) as f(s) in the 
half-plane > ^ . This completes the proof of Proposition [TJ □ 

Remark 2. As a direct application of the theory developed in Jessen [12] and |13j . 
we also obtain the following facts regarding the zeros of any Dirichlet series f(s) 
satisfying all assumptions of Lemma [TJ Let f(cr) be the Jensen function corre- 
sponding to /(s) as introduced in [12]. It follows from [121 Satz A and p. 492] and 
[T3l §28] (applied to the equation ]T°1 2 e ( Al " A ^ s = -1) that (p(a) is a C 2 -function 
with non-negative second derivative throughout K >(J0 . For any fixed o\ < 02 in 
^>o- ) we let A r (o"i, 02, T\, T2) denote the number of zeros of f(s) in the rectangle 
s G (<7i,0"2) x (T\,T2) (counting multiplicity); then 

(2-4) hm = — / ip (a) da. 

Ti-T^co I2-I1 27T J ai 

Furthermore, it follows from [12, Satz B] and Lemma [TJ above that (p"(cr) does not 
vanish identically on any subinterval of (00,07), ^ ne nm it m ([231) i s positive for 
any fixed o\ < 02 in M>o- with o\ < Of. 

In particular this applies, via (|1.4p and Lemma [21 to the zeros of E n (L,s) for 
almost any L € X n . Thus, if L is generic in the sense of Lemma [21 then we conclude 
that for any fixed o"i,o"2 with § < <ti < 02, the limit relation (|2.4p holds with 
N(ai, CT2, Ti, T2) now denoting the number of zeros of E n (L,s) in s 6 (o"i,<T2) x 
(Ti,T2); furthermore the limit in (|2.4p is positive whenever a\ < ai). 

In a recent paper by Lee [16], similar asymptotics for the number of zeros are 
obtained in the more difficult case of a rational L G X2 corresponding to an integer 
quadratic form with a fundamental discriminant. 

Remark 3. Consider the function 



L i->- a L 



a > § such that 2m(Ly 2(T = \C{2o-y l E n (L , a) 



(Thus Proposition [JJ says that ul = almost everywhere.) Let X' n be the (closed) 
subset of X n consisting of those lattices for which # {v G L : \v\ = m(L)} > 2; by 
|28l Lemma 5.1], X' n has measure zero. We claim that ox is a smooth function from 

X n \ X' n tO 

^>n/2- This follows by studying the following function on K >n / 2 X X n : 
a(a,L) = 2m{L)- 2a - ^C(2a)- 1 E n (L, a) = 2m(L)- 2a - ^ \v\- 2a . 

This function is smooth on all M> n /2 x (X n \X' n ) and one easily checks that -^a(a, L) > 
at all points a = &l, L G X n \ X' n . Hence our smoothness claim follows from the 
implicit function theorem. 

On the other hand note that ai — > 00 whenever L — > Lq for some Lq G X' n . In 
particular this shows, via Proposition [U that sup Lg ^ n <tl = 00. (But of course, as 
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we remarked in the introduction, 07, is finite for any fixed L G X n , in particular for 
any L G X' n \) 

Remark 4. Note that for any rational L G X n , there occur arbitrarily large multi- 
plicities among the lengths \v\ for v G L; in particular, the independence condition 
in Lemma [2] fails for every rational L G X n . Indeed, for n > 3 this claim follows 
easily from the fact that the number of v G L with \v\ < R grows like R n as R — > 00, 
while the number of possible values of |«| grows at most like i? 2 (since L rational 
implies that there is some c > such that |i>| 2 € cZ for all i; G L). The claim also 
holds for n = 2, since in this case the number of possible values of \v\ with \v\ < R 
is in fact < Pt 2 (logPt)"5 as i? -> 00; cf. [2] or [18]. 

3. Proof of Theorem [TJ 
The sequence {Tj}'jL 1 of points of the Poisson process V belongs to the space 



n 



x = { Xj }f =l G (R >0 r : < x x < x 2 < x 3 < . . . }, 



which we equip with the subspace topology induced from the product topology on 
(M>o) °. We denote the distribution of V on ft by P; this is a Borel probability 
measure on ft. 

Recall the definition (jl.6p of 0~{t }> let us prove that this is a well-defined random 
variable on ($7,P). We set 



x = {xj}f =1 G ft : #{ Xj < X} ~ ±X as X -»• 00}. 



This is a Borel subset of ft and, by the strong law of large numbers, we have PC = 1. 
For any x G O', we have Y^jLi x J 2a < 00 f°r an > 5 ana ^ S^Li x 7 2fT ~~ ^ 00 as 
a — > and thus, by the same argument as in the proof of Lemma [Tj there exists 
a unique a > \ satisfying x\ 2(J = Y^f=2 x ^ 2<T ■ ^ n other words, o x is well-defined for 
every x G ft'. Furthermore, for any c > |, we have {a; G f2' : a x < c} = {x G ft' : 
x^ 2c — X^2 x 7 2c > 0}, which is a Borel set. This proves that the function x 1— > a x 
is P-measurable on ft, i.e. that o"{t } is indeed a well-defined random variable. 
For given n > 2 and c > ^, we let F n (L, c) be the random variable given by 

00 

F n (L, c ) :=-Vi(i:r 2c + X>(i:r 2c , 

where as usual L is taken at random in X n according to fi n . We also let F(c) be the 
random variable 

00 

F(c) := -Tf 2c + ^rr 2c . 

i=2 

Lemma 3. Lei c> \ be fixed. Then F n (L,c) converges in distribution to F(c) as 
n — > 00. 

Proof. The proof is a straightforward adaptation of the proof of [27, Thm. 1] with 
m = 1. □ 

Lemma 4. For any given c > \ and r£l, we have P{-F(c) = r} = 0. 
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Proof. The lemma follows immediately from the calculations in the second paragraph 
of the proof of Theorem [2] below. □ 

Proof of TheoremUl It suffices to prove that for any fixed c > ^, Prob^ n {n^ l aL > c) 
tends to P(<7{^.} > c) as n — > oo. By Proposition [1] and the monotonicity argument 
at the beginning of the proof of Lemma [H we have 

oo 

ProV>-V L > c) = Prob^L € X n : -2Vi(L)" 2c + C(2cn)- 1 ^ V 3 (L)- 2c > o} 

oo 

= Prob^L € X n : F n (L,c) > (l - C(2m)~ 1 ) ^ V^L)-^}. 

i=i 

Now let e > be given. By LemmaHwe have P{-F(c) = 0} = 0, and hence (using 
[191 Thm. 1.19(e)]) there exists r > such that 

(3.1) P{F(c) e [0,r]} <e. 

Furthermore, it follows from |27| that there exists K > and iV £ Z + such that 



oo 

Pro6^{L € X n : ^ V,(L)- 2c < > 1 - e, 



Vn > AT. 

i=i 

After possibly increasing N, we may also assume that (1 — Q{2cn)~ 1 )K < r for all 
n > N. It follows that, for all n > N, 

Prob^ n {F n {L,c) >r)-e< Prob^n^aL > c) < Prob^ n (F n {L,c) > 0). 

However, by Lemma [3] and Lemma HJ we have (cf. [3l Thm. 2.1(v)]) 

lim Prob^ n (F n (L,c) > 0) = P(F(c) > 0) = P(^ {T?} > c) 

n— >oo v 

and 

lim Prob^ n {F n {L,c) > r) = P(F(c) > r). 

n— s>oo x 7 ' 

Furthermore, by ([3TT]) we have P(F(c) > r) > P(F(c) > 0) - e = P(cT{ T } > c) - e. 
Hence we obtain 

lim sup Prob lln (n~ 1 o- L > c) < P((T{t j } > c ) 

n— >oo 

and 

lim inf Probn n (n^ai > c) > P(err T .i > c) — 2e. 

n— >oo y v \ JJ / 

But e is arbitrary and hence the proof of Theorem Q] is complete. □ 

4. Proof of Theorem [2] 

In this section we prove Theorem [21 which gives an explicit formula for the distri- 
bution function of 0"{t.,}- Recall that cr^} > c holds if and only if Y^=2 Tj~ 2 ° > c . 
Hence our task is to determine the probability 

oo 

(4.i) p(^ } > c )=p(E :z 7 2c > r r 2c )- 

3=2 

Let us note that for any fixed /i > 0, the sequence fj,T±, //T2, . . . give the points 
of a Poisson process on the positive real line with intensity (2fi)~ 1 , and we have 
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S^=2^j~ C > c if and only if X^UO^j) -20 > (A^l) -20 - Hence we may, in order 
to make our computations slightly cleaner, alter our notation so that from now on, 
< T\ < T2 < . . . denote the points of a Poisson process on the positive real line with 
constant intensity one; the probability in (|4.ip remains unchanged by this alteration. 
Also to make the computations slightly cleaner, we will write 

a := 2c € R>i. 

As a first step, we consider the conditional distribution of the sum ^^=2 Tj~ a 
given the value of T±. We will see that this distribution is infinitely divisible. For 
basic facts about infinitely divisible distributions, cf., e.g., [TDJ, Chs. VI. 3, IX, XVII]. 
We formulate the result for a Poisson process having constant intensity 1; it is of 
course easy to carry this over to the case of an arbitrary constant intensity. 

Proposition 2. Let < T\ < T2 < ■ ■ ■ be the points of a Poisson process on the 
positive real line with constant intensity 1. Then, for any a > 1 and 5 > 0, the 
conditional distribution o/^^L 2 ^~ a > given that T\ = 5, is an infinitely divisible 
distribution, the characteristic function of which is given by 

(4.2) <p a>s (t) = E^e u ^ T r I Tx = 5) = expj- jH(l - e Ux ~ a ) dx\. 

(Cf. [20, Thm. 1.4.2], where the corresponding fact is proved in the special case 
5 = but with more general weights in the sum; the resulting distribution is then a 
stable distribution.) 

Proof. Let n be a positive integer and let r\ be any real number larger than 5. 
The conditional distribution of (T2, . . . ,T n+ i), given that T\ = 5 and T n+ 2 = r], is 
that of the order statistic of n i.i.d. random variables uniformly distributed in the 
interval (5, rj), and hence the conditional distribution of Y^=2 Tf a \ given T% = 6 and 
T n +2 = T], is the same as the distribution of Y^j=x(^ + (v ~ &)Uj)~ a -, where from now 
on Ui, U2, ■ ■ ■ denotes a sequence of i.i.d. random variables uniformly distributed in 
(0,1). It follows that the conditional distribution of Sj=2 -^V"' gi ven onr y ^1 = ^ 
is the same as the distribution of 

n 

X n :=Y J (S + S n+ iU j )- a , 
3=1 

where S^+i denotes the sum of n + 1 i.i.d. exponential random variables with mean 
one, independent from the sequence {Uj} (so that 5 n+ i has the same distribution 
as T n+2 - 5 given Ti = 5). 

By the law of large numbers n S n +% tends in distribution to 1, i.e. given any e > 
there is £ Z>o such that for each n > N, we have (1 — e)n < S n +i < (1 + e)n 
with probability > 1 — e. It follows that if we let 

n 

Y n :=Y,{5 + nUj)- a , 
i=i 

then, for each n > N, we have (1 + e)~ a Y n < X n < (1 — e)~ a Y n with probability 
> 1 — e. In particular, it now suffices to prove that Y n tends in distribution to a 
(non-defective) random variable whose characteristic function is given by the right 
hand side of f|4.2[) . since then also X n must converge in distribution to this random 
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variable, and also it follows from the definition of Y n that the limit distribution must 
be infinitely divisible, cf., e.g., [TDJ Ch IX.5 (see also Ch. XVII.2)]. 

But Y n is a sum of n independent random variables, and thus its characteristic 
function equals 

Ee UY n = ^(S+nU^y = A J* +U e to- d X = Ll J^ +n ^ _ git*-) dx 

Note that |1 — e ltx a \ -C |t|a;~ a uniformly for all x > 5 and all t G E. In particular, 
for each fixed t £ I the integral J/°(l — e ttx a ) dx is absolutely convergent, and 
Ke itYn tends to the expression in the right hand side of (|4.2p as n — > oo. The bound 
|1 — e ltx a \ <C |i|x _a also implies that the function (p a ,s{t) is continuous. Hence Y n 
converges in distribution to a (non-defective) random variable whose characteristic 
function is given by the right hand side of (|4.2f) . and the proposition is proved. □ 

Remark 5. Let us note that the integral in (|4.2p may be expressed in terms of the 
incomplete gamma function. Indeed, substituting x = (iit)" and then integrating 
by parts, we get 



oo 



l-e^ a )dx = -r (l-e-»)(±(W- k *))du 



it&~ a -i\ i . / „-tur - 1 



= 5{e lW -l)+t / e- tu {iu)-a du 
Jo 

(4.3) = 5(e its ~ a - 1) + (-i*)- 7 (l - \, ~it5~ a 

where for t ^ we agree that arg(— it) = — (sgnt)^. Hence 

Va , 5 (t) = exp{-<5(e^~ a - 1) - (-it)* 7 (l " "^" a )}- 

Furthermore, using the recursion formula (|1.8p together with the formula j(s, z) = 
T(s) — T(s,z), where 

/oo 

is the upper incomplete gamma function, we get the alternative formula 
(4.5) <p a , s (t) = exp{5 - (-it)lr(l - i) - i(-zi)ir(-i, -rt^)}. 

Proof o/ Theorem^ Note that, for all z <E C \ {0} with J?z > 0, 

|r(-i,*)|<w-i- 1 e-«'. 

Hence, if we denote the exponent in (|4.5p by ipa,s(t), we have for i > 0, 

^M*) = ^ - H*)*r(i - i) + o M (r 1 ) . 

Using also it)« 3> a we conclude that —$top a ,5{t) ^a,8 t° as t — > oo. Hence, 
in view of the symmetry (p a ,s(—t) = Pa,s(t), the function (p aj s is integrable, and 
therefore the distribution in Proposition [2] has a density function, which we call 
fa,s(x). Thus 

-i roo i roc 

fa, S ( X ) = n - / Va,8(t)e~ itX dt=- ^{^ S {t)e- itx ) dt. 

^ J-oo ^ Jo 
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It follows that the conditional probability of J^^-^T* > T^ a , given that T\ = 5, is 

OO n QQ 

P(VTp > Tf a \T 1 = S)= fa, s (x)dx. 

K U 1 Js - a 

However, Ti, being the first point of a Poisson process on the positive real line with 
intensity one, has an exponential distribution of mean one. Hence we conclude: 



OO n QQ /"OO 

{a {Tj} > c ) = p(E T r > T r a ) = / /_ dxe -~ s dS 

o 

! f a ,s{x)e~ 5 d5dx 



3=2 



IT 



X a 

oo />oo 



(4.6) =-/ / / R(cp a>5 (t)e- itx - s )dtd5dx 



X a. JO 



Note that the last expression in (|4.6p should be viewed as an iterated integral; it 
is easy to see that Jq 00 J°1i / °° \$t{<Pa,s(t)e~ ltx ~ s ) | dtdS dx = oo, so that we are not 
permitted to change order of integration arbitrarily. However, we will prove that the 
inner double integral is absolutely convergent. 

By Proposition [2] we have e~ s ^p a ,&{t) = exp(— (5 + — e ltx a )dx)), and here 

we have, by substituting x = (u/t)~° and then integrating by parts, 

5 + J {l-e itx ' a )dx = 5 -t^ J (l-e iu )^(u-^du 

= 5e it6 ~ a - it* / e iu u-^ du = t^ a (t5~ a ), 
Jo 

where we have defined 

i*y 

(4.7) ® a (y) := y~*e iy - i e %u vT^ du for a > 1, y > 0. 

Jo 

Thus 

1 /-OO /-OO /-DC 

(4.8) P(a {T . } > c) = - / / / Kexp|-ite-^$ a (t,5- a )|dtd<5dx. 

Using e iu = 1 + O(u) for u G [0,1], we find that $ a (y) = + O(y)) for 

< y < 1. (Here, and in any "big-O" or "<C" bound below, we allow the implied 
constant to depend on a.) In particular there exists a positive number K\, which 
may depend on a, such that 5M> a (y) > \y~~^ for all y G (0, ki]. We also note that 

(4.9) K(y) = - -y- 1 -*^. 

a 

In particular K<3?„(y) = — \y~ X ~ » cos y, and this is negative for all y G (0,^), so 
that 3?<I> a (y) > K^ad) holds for all y G (0, §). Furthermore, for all y > f we 

have 9?$ a (y) = K# a (f ) - ± J| u -1 ~a(cosu) du > K$ a (§). Also note from gTTD that 

2L 1 

5? < I > a (|) = Jq 2 (sinu) > 0. Hence we conclude: 

$t$a(y) > *2 := K* a (f ) > 0, Vy > 0. 
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e~s 5 dt + / e~ K2ta dt < e~ Kz 

J Kx5 a 



Using the bounds obtained, we conclude: 

(4.10) J exp{-^$ a (i<T a )} dt < J 

for all 5 > 0, where K3 is some positive number which may depend on a. From this 
estimate we see that the inner double integral in (|4.8p is indeed absolutely convergent, 
in fact even / °° J °°|exp(— ta<& a (t5~ a ))\ dtd5 < 00. Hence we have 

1 l-X t-OD t-OD j 

P(cr {Tj} > c) = - lim K / / / x ex.pl -itx -ta^ a (t5~ a )\ d5dtdx 



- lim 3? 

7T X— »oo 



00 ,>A 



JO 

00 , 



= — lim 5? / t* 

ira x-5-oo Jq 

Here, for any t > 0, we have, by integration by parts: 



-itx 



X 



[ 1 expj-^$ a (t(T a ) j dddxdt 

J X~ ~ 

e -t-^(y) y -l-l dy dx dt 



x 



-itx 



ie 



-itX rtX 



X 



tx 



-it Xe -d^ a (tx)^ x yl-\ dx 

<i/) e -*«*.(w) y -i-i 



Hence 



P(a {T . } >c) 



— lim 9 / t*" 1 / (e-^-e-^e-^^^y-^sdydt. 



For given X > 1, we split the integral over t into two parts, corresponding to 
t < X^ 1 and t > X^ 1 . Regarding the first part, we note that t < X~ x and y < tX 
implies y < 1. Thus $> a (y) = J/~"(l + 0(y)), and since also t < X -1 < 1, we 
have e"*"*' 1 ^ = e" ( * /s/)7I (l + O^y 1 -^)). Recall that in this case we also have 



-ni 



1 + 0{y). Hence 



(4.11) 



x- 



--1 



tx 



e~ iy -e 



x- 



IX 



uxy-tlt*a<3i)y-i-± 
l-e~ itx + o(y + tay 1 -* 



dy dt 



•(*/v)«^-i„-i-4 



t" y a dy dt 



X' 



i (1 - e~ ltx + 0(tu" a + ^ 1 - a ))e~ u t- 1 dudt 



1 j _ e -« 



dt 



1 

X~a 



e- u du + 0(X 



1 

X~a 



[ u - a + u 1 - a ) e - u du 



l-e 



The remaining part is 



dt + 0(X~). 



(4.12) 



--1 



A- 



tX 



dy dt, 



14 ANDREAS STROMBERGSSON AND ANDERS SODERGREN 

and here we have absolute convergence; Jx-i Jq X e~' a5R<I>£l( ' ? '^y _1 ~a dy dt < oo, 
as is seen by a similar computation as in (|4.10p . (The corresponding fact does not 
hold in (|4.1ip .) We also note that we may replace the range of the inner integral in 
(|4.12|) by all of R>o, to the cost of an error which is 

/•oo />oo i />oo i if 

</ ta- 1 €-^ ta y- x —adydt^X— a \ e~ K2ta - < X~ log(2X). 
Jx- 1 Jtx Jx- 1 t 

Collecting the above results, and using the fact that both X~a and X~a log(2X) 
tend to zero as X — >■ oo, we conclude that 

(4.13) 

p ('m)>«) = ;/^* + i J fi5.(C s ««W*-jr i s »(«- 1 "»'W)* 



where 

9£ 



i r°° i 1 

(t)=t^- 1 / e-^-^My- 1 -* dy 
Jo 

for £ = 0, 1. 

Next, we split g e (t) as g e (t) = ge,i(t) +gt,2{t), where 

g lA (t) = t 1 *- 1 ! l e-««-***-(v) y -i-i dy 
■/ o 

and 



Bounding 3?$ a (y) from below as in (|4.10p . we see that for all t > we have 
(4.14) < i^ 1 / 1 |e~* i ' My) |?r 1 ~" rfl/ < rV 4 *", 



o 



where K4 is (just like k\, K2, K3) a positive number which may depend on a, and 

/OO 1 1 



Note also that for all t,y G (0,1], we have e - * 1 ^* *-^ = e~ iy - tay a ( 1 +°^ 

e -tay~a ^ + ^I y l-I + ^ ^ thus 

ia -1 /" I Qfe - * 1 ^-*" [ J/~ 1_ » dl/ < t^ 1 f (fry- 1 * +y~l)e~^ y ^ dy 
Jo Jo 

/•oo 

(4.16) < / (v 1 - 11 + v- a )e~ v dv ^.T*- 1 

Jta 

for all < t < 1. Combining this bound with (|4.14|) and (|4.15p . we see that 

/■oo />oo ^ 

(4.17) / t~ l ^e-^-^^ly- 1 - a dydt< 00. 
Jo Jo 
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Hence the contribution from gi(t) in (|4.13p can be treated as follows: 
(4.18) 

1 f°° 1 f°° i f°° i i 

— lim / $ gi (t)dt = — t" 1 ^e- iy - ta ^ y) y- l -«dydt 
ira x^oo J x -i ira J J 

Finally, we treat the contribution from go(t) in (|4.13p . Note that, by (|4.14p and 
(|4.15p . the restriction of go(t) to [l,oo) is an L 1 -function. Hence, by the Riemann- 
Lebesgue lemma, e~ tXt go(t) dt tends to as X — > oo. Moreover, the restriction 

of 5o,2(*) to (0, 1] is in L 1 and hence also f x -i e~ lXt go,2(t) dt tends to as X — > 00. 
Hence 

1 f'°° 1 f' 1 

(4.19) lim / %(e- lXt g (t)) dt = lim 9 / e~ lXt g 01 {t) dt. 

7ra x^oo J ' 7ra A->oo Jx- 1 

Furthermore, for < t < 1, we have 



a 
1 



where we bounded the contribution from the big-O-term in the integral by a similar 
computation as in (|4.16p . Thus go,i(t) — f is an L 1 -function on t 6 (0, 1], so that 
J^_i e~ iXt (g ,i(t) - j ) dt tends to as X ->■ 00. Hence (|4Tl9|) equals 

1 . „ f 1 e~ iXt , 1 . f x sint , 1 f^sint 



dt = - lim / dt = - dt. 



7T X— s-oo Jx -1 t ft X^qci Jit IT Ji t 

Collecting our results into ()4.13p . we obtain, since J °° dt = ^, 

< 42 °> p ('™>' : )-5 + ?r 9 (^) iri " i *- 

Let us note that <& a (y) can be expressed in terms of the incomplete gamma function, 
by substituting u = iv in (|4,7p and using formulas (|1.7p and (jl.8p : 

1 1 — j-i 

(4.21) $ (y) = y^e*" + e^Vfl - -, -iy) = -^-^yf — > ~% 

V a / a \ a 

Substituting this into (|4.20p , we obtain the formula stated in Theorem [2j Using 
|$«(y)| > K$ a (y) > k 2 > for all y > and $„(y) = y~i(l + 0(y)) for < y < 1, 
one immediately sees that the integral in (|4.20p is absolutely convergent (this is also 
clear from the proof, cf. in particular (|4.17p and (|4.18p ). This concludes the proof 
of Theorem [2J □ 

Remark 6. It is worth stressing that if we remove the imaginary part in (|4.2U|) . then 
convergence fails: We have \J yo £ ¥^ y~ l ~a dy\ — > 00 as yo — > + , since 3> a (y) = 

y~u(l + 0{y)) for < y < 1. 
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5. Proof of Corollary Q] 

In this section we prove Corollary [TJ To begin, note that by formal differentiation 
under the integral sign in (|4.20p . we have Prob[a^x 3 } < c) = Jy 2 /(ci) dc\, where 
/ : IL i — > M>o is given by 

Here a := 2c G R>i (see Section This manipulation is justified by the fact 
that the integrand in (|5,ip is majorized, uniformly for a in compact subsets of K>i, 
by an integrable function; this follows from an argument similar to the one that 
shows that the integral in (|4.2U|) is absolutely convergent, using also that -§^&a(y) = 

a~ 2 {\ogy)y-- a (1 + 0(y)) for < y < \ and £$ a (y) = O(l) for \ < y < oo. 

Remark 7. Note in particular that the imaginary part in (|5.ip may be taken outside 

the integral; in fact even / °° I I^J ( %^(yf ~ l? 2 ) \ y dy < 00 ■ 

Let us now consider formula (|5.1|) in the limit as a — > oo. In (|4.7p . we expand 
e tu in a power series, change order between summation and integration and then use 
(n — a -1 )^ 1 = n^ 1 YlT=o( na )~ k f° r eacn n G This gives 

oo 

(5-2) <S> a (y) = y-- a (l-Y J F k (y)a- k 

k=l 

where 

" iiy) n 



(5-3) F fc (y) := ]T 



n=l 



Obviously |i*fc(y)| < e' y ' — 1 holds for all y > and all k, and hence we see that given 
any yo > there exists some ao = ao(yo) > 1 such that £)fcLi Fk{y)( 1 ~ k \ < 5 holds 
for all a > ao, y € (0,yo]. We also have |X)fcLi Fk{y)a~ k \ *C |y|a _1 for these a,y, 
and therefore $ a (y) _1 = y«(l + F\(y)a~ l + 0(\y\a~ 2 )). The power series in (|5.2p 
may also be differentiated termwise with respect to a. Using these observations, we 
obtain by a short calculation: 

(5.4) VC^ f tMy) \ogy\_F 1 (y) | 2{F 1 (yf + F 2 (y)) | Q fl + \logy\ 



V *a(l/) « 2 / ya 2 ya 3 V a 4 

uniformly over all a > ao(yo), y £ (0, yo] (where we recall that yo > is arbitrary). 
In order to obtain a similar relation also for large y, we start by setting 

(5.5) e a := lim <D a (y) = ^ 

2/->oo a 

(cf. (jOTj) ). In view of we have $ (y) = f(a) + a" 1 J^ 00 u~ 1 ~^e iu du, and 

integrating by parts twice, we get (for any a > 1, y > 0) 

(5.6) $ a (y) = £(a) + ^r(0, -iy) - ^ll(y) + -I /°° u" 1 -^^) (in, 
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where 

,5.7) n ,, ); =rn^) d „. 

Jy U 

We have |r(0, — iy)\ <C y~ l for all y > 0, and thus also |n(y)| <C y _1 for y > 1. Using 
this fact together with the trivial observation — 1 — ~<— 1, we bound the integral 
in (|5.6p and get 

_i _i 

(5.8) d> a (y) = £(a) + ^r(0, -iy) - ^n(y) + O^y- 1 ), 

a a z 

uniformly over all a > 1, y > 1. Since also £(a) = 1 + (7 — ^a" 1 + 0(a -2 ) as 
a — > 00, we see that <& a (y)/£(a) is near 1 whenever a and y are large; hence there 
exist absolute constants a\ > 1 and yo > 1 such that for all a > a\ and y > yo, 

_ j_ 

' ' ° r r(0, -iy) + 0(a-V 1 ) 



2" 



<My) £(a) a£(a) 

(5.9) = ^_^r(0,-iy) + O(a-V 1 )- 

£(a) a 

In order to obtain an asymptotic formula also for J^$ a (y), we note that the right 
hand side of (|5.6p defines an analytic function of the complex variable w = a -1 in 
the region \w\ < 1 (including w = 0). Restricting to \w\ < |, we may bound the 
absolute value of the integral in (|5.6p using [n(u)| <C u- 1 and K(-l-w) < -\. We 
may then use the Cauchy differentiation formula to obtain an asymptotic formula 
for the derivative of our analytic function, valid uniformly for \w\ < \. In particular, 

(5.10) 

— $ a (y) = £'(a) - ^r(0, -iy) + ^(W(y) + T(0, -iy) logy) + 0[ar*y-*), 

uniformly over all a > 4 and all y > 1. Using (|5.9p and (|5.10j) . we obtain, via a 
straightforward computation, 

(5.H) i 

y'^fm^aiv) logy\ _ £» _ (log y)y _1 -- y~ 1_ «r(0, -iy) 



$ a (y) V *a(l/) « 2 / £(a) 2 « 2 ^(«) « 2 

+^^|(4 7 - 2vri + 21ogy + 2y"^r(0, -iy))T(0,iy) + 2n(y)| + 0( a - 4 y-§), 

uniformly over a > max(ai,4) and y > yo- 

We now multiply the relation (|5.1ip with e~ iy , and integrate the result over y € 
[yo, 00). The contribution from the first term is 

(5.12) g| r v- 1 -**-* *y- 



We split this integral into two parts as fy* P ^ a 1 ^ + /exp^a 1 /*) (keeping a so large that 
exp(a 1//4 ) > yo); then, because of the oscillating character of the integrand, the 
second integral is 0(exp(— a 1 / 4 )). In the first integral, we use y~« = 1 — + 
i(l2p)2 + ((l2M)3) a nd (iy < a; then, by a quick computation, we 
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find that ([Ell equals fg$(/~ dy- \ J™ 2l!^2£l dy + 0( a - 2 )). The remaining 
terms in (|5.1ip can be treated similarly, and using the relations 

(5.13) Fi(y) = ft- 7 -logy-r(0,-iy); 

F 2 (y) = II(y) + (^vr 2 - ± 7 2 + ±7177) + {\m - 7 - | logy) logy, 

the result may be collected as 



2 



=a -2 rEM e -iv iy+a -3 r w<*r+w)f+ iy+0{a - % 

J ya y Jyo y 

for all a > max(ai,4). Using also (|5.ip and (|5.4p . we thus obtain an asymp- 
totic formula for /(c) as c = \a — > 00. Note, however, that G J °° e~ ty dy = 

\^ I^oo Fl y V ^ e ~ iy dy = 0, where the second equality follows using the Cauchy inte- 
gral theorem, moving the contour towards infinity in the lower half-plane. Hence 
the coefficient in front of a -2 = (2c) -2 in the asymptotic formula vanishes, and we 
arrive at ()1.10p . with 

(5.14) K 2 = —9 r Fl ^ 2 + F2 ^ e -iy dy = 0.822467 .... 

2?r J y 

(The numerical evaluation of this integral, which is not entirely straightforward, is 
carried out in [25] constants. mpl].) 

We next turn to the study of (|5.ip in the limit as a — > 1. Our presentation here 
will be rather brief; we refer to |25^ asymptotics.mpl] for further details. The formula 
(|4,7p may be expressed as 

(5.15) $ a (y) = ySjy - ie iy -^j - — ^ f e iu u l -~« du. 

1 - a 1 1 - a 1 In 



Now fix N G Z+, and let us keep (a - 1) N < y < (a - 1) -Af , and a G (1,2]. We 
split the integral in (|5.15p as /J a ^ + JZ^jv+i and bound the first part trivially, 

while for u G [(a — l) N+1 ,y], we use the fact that u 1 ^^ = ^^Jq 1 ^— - — fcp^^ 
+0((a — l)^ 2 ! logu\ N+2 ), where the error is an increasing function of u when u > 1. 
This leads to the formula 

(5.16) * a („) = r -^ T 1 + £ G»(y)(l - O fc + o((a - (l + -)) , 

^ k=i y ) 

where G\(y), ^(y), . . . are given by 

Let us now further restrict to the case where (o — 1)2 < y < (a — 1) -Ar . Using 
\Gk{y)\ I l°g y| fc X y 1 + I l°gy| fe + 1, we see that there is some ao = ao(iV) G (1,2] 
such that for all a G (l,ao] and all y G [(a — 1)2, (a — l) _Ar ] the expression within 
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the brackets in (|5,16p lies in {z : \z — 1| < i}, and so we get 
1 f N ( - 1 £ 

— .^-Oji+Ej-XW-.- 1 )*} 

+0(( 9 - 1 + |log 9 |) A ' +1 (a-l)"+')}. 

Working similarly, starting from a differentiated version of (|5. 15[) . we also get an 
asymptotic formula for J^$ a (y), an d with further computation, we finally obtain 

-l-i / 8 f..\ !„.,„. \ ■ -i V ( N 



My) V *a(y) a 2 / y 2 I ^ 

+o((y- 1 + |iog y |) Af+1 (a-i) iV+1 )}, 

for all a € (l,ao] and y € [(a — 1)2, (a — 1) _JV ]. Here Hi(y), H%(y), . . . are certain 
continuous functions of y satisfying |i?|(y) | <Q (y _1 + | log y\) e ; in particular we have 

(5.17) 

H 1 (y) = 2U J (\ ogu y u du- — + (l- eulogy 

( i f y n ■ ( f y i e iy ■ 

H 2 (y) = 3 1 - J (log u) V" du - I J (log u)e iu du + i - - log y - — + ie JJ/ log y 



ie iy log y 



logy+(i-ie^)(logy) 2 |. 



y 

Writing Ht{y) := —ie~ iy y~ 2 Hi(y), it follows that, for y < 1, 

Qi^y) = 2y~ 2 + i - ^y 2 + 0(y 4 ); 3ff 2 (y) = 3y" 4 + fy" 2 - ±? + 0(y 2 ). 
Furthermore, one computes (again for y < 1) 

G# 3 (y) = -4y" 4 - f y- 2 + 0(1); 9F 4 (y) = -5y' 6 - fy" 4 + 0(y~ 2 ); 

ZH 5 (y) = 6y- 6 + 0(y~ 4 ). ZH 6 (y) = 7y" 8 + 0(y" 6 ). 

Using these relations (taking N = 6), we obtain 

J(a-i)i *a(y) V *.(y) « 2 r y "io y 2 y 

+{a - 1) y (^1 (y) - 2y- 2 ) dy + (a - l) 2 jf (0fH 2 (y) - 3y~ 4 - |y~ 2 ) dy 

(5.18) +{-(a - I)"! + |( a - 1)1 - f|(a - 1)1 - ™(a - 1)* } + 0((a - l) 3 ). 

(This formula is first derived with each upper integration limit being (a — l) -4 (say) 
in place of 00; the remaining integrals over y £ [(a — 1)~ 4 , 00) are easily seen to be 
subsumed in the error term.) 

To treat the integral over y < (a — 1)2, we start with the formula 
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which holds uniformly over all a > 1 and < y < 1, for any fixed N € Z>2; this is 
proved using (|4.7p and the power series expansion of e m . Note that the sum over k 
is 0(y min(a — 1, y)); hence there is an absolute constant yo £ (0, 1] such that for all 
a > 1 and < y < yo, we have 

l/°(o-l) L , /y^ i fc (a-l) 



*o(y) a-l-iy 



\<i<N/2 v fc=2 



Using this formula with N = 5, together with a similar asymptotic formula for 
~§g^a{y) deduced from a differentiated version of (|4.7|) . we find after some computa- 
tion that 

e~^~ 1- ° (-esPM logyA _ P (a-l ) y) + Pi(a-l,y)logy 



*a(y) V *a(l/) a 2 / a 2 (2a-l) 6 (3a-l) 5 (4a-l) 4 (5a-l) 4 (a-l-%) 

+o((a-l)y 3 (l + (a-l)|logy|)^ 

where Po an d -Pi are explicit polynomials. This formula can now be integrated over 
y in terms of elementary functions, and we obtain 

Jo *a{y) V « 2 r y " 27r 2 " (a l) 

(5.19) +{(a - - f (a - 1)4 + f§(a - 1)1 + ™(a - 1)*} + 0((a - l) 3 ). 

Finally, we add (ISTTHj) and (15331) . and note that since Pi(y) = -2(Pi(y) + ^ + 1), 
we have 

J°° [QHx (y) -2y~ 2 )dy= [°° ^ ^fi^ dy - vr = 0, 

where the second equality follows by again moving the contour towards infinity in 
the lower half-plane, noticing the pole at y = 0. Hence we arrive at (|1.9[) . with 

(5.20) ^i = 20 + -y (%~ 4 + §y~ 2 -%H 2 (y) J dy = 39.47841... 

(cf. [25^ constants.mpl]). This completes the proof of Corollary [TJ □ 

Remark 8. It appears that by the same method one could obtain asymptotic expan- 
sions of /(c), in the limits as c — > 00 and c — >• i, with the error term having any 
desired power rate of decay. 



Appendix A. Residue calculus and numerical computation 

OF THE DENSITY 

In this appendix we discuss the evaluation of the integrals in (|4,20p and (|5.ip using 
the residue theorem, resulting in alternative formulas for P(cj|jv} > c) and the cor- 
responding density. These formulas turn out to be useful for numerical computation, 
something which we discuss briefly towards the end of the appendix (see also [25] 
numdensity.mpl] ) . 
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We now write z in place of y. By (|4.7p we have & a ( z ) = z a {e tz — iz e lzt t « dt), 
and here the expression in the parenthesis is clearly an entire function of z. Hence 

(A 1) * (z) — e ~ lZz ~ 1 ~^ — e~ %z z~ l 

$a{z) e iz - iz Jq e izt t~a dt 

is a meromorphic function in all of C. In (|4.20p we are integrating Q^^z) along the 
positive real line; using the symmetry 1 3/ a (— z) = —^ a (z), we may rewrite this as 

(A.2) P(a {Tj} >c ) = \+ito±( 12 *.(y) dy + f Mv) dy) . 

Let C' r be the semicircle {z : \z\ = r, Qz < 0}, oriented in the direction from —r 
to r, and let C r be the contour going from — oo to —r along K, then from — r to 
r along C' r and finally from r to +oo along R. Since ^ a ( z ) has a simple pole at 
z = with residue 1, we have f c , ^ a (z) dz = m + 0(r) as r — > 0. Thus (|A.2|) equals 

lim r _ >0 + 2^? Jo x ^a( z )^ 2: - However, by Cauchy's integral theorem, J c ^ a (z)dz is 
independent of r for all sufficiently small r. Hence 

(A.3) P(a {T . } >C )=^-J c * a (z) dz, 

for any r > so small that ^ a {z) has no pole in the punctured disk {z : < \z\ < r}. 

We wish to replace C r in (|A.3[> by a contour over z's with large negative imaginary 
part. In order to do so, we first need to understand the poles of *& a ( z ) m the lower 
half plane. Numerics indicate that there is exactly one simple pole in the infinite 
vertical strip {z : (2n — l)n < $tz < (2n + 1)tt, Qz < 0} for each integer n; cf. Figure 
[2] below. However, for technical reasons it seems easier to prove a corresponding 
statement instead for certain "curved vertical strips", as follows. For each n € Z + , 
we let T n be the curve in the complex plane given by 

(A. 4) x i->- Cnix) = x — ixt&nUn — j)tv — hx) , (2n — |)7r < x < (2n — 

One notes that ^sc n (x) — > — oo as x — >• (2n — |)7r + , that Sc n ((2n — |)7r) = and that 
< argc4(x) < | for all (2n — |)7r < x < (2n — h)ir. Hence T n and r n+ i, together 
with the real interval [(2n — ^)ir, (2n + |)vr], bound a curved vertical strip, which we 
call S n (we take S n to be closed). We also let S- n = {—~z : z £ S n } be the reflection 
of S n in the imaginary axis, and we let So be the curved vertical strip bounded by 
the curves Ti, {— z : z G Ti} and [— fvr, |-7r]. Now the union of all S n (n € Z) equals 
the negative half plane, {z : Qz < 0}, and the SVi's have pairwise disjoint interiors. 

Proposition 3. Let a > 1 be given. For each n E Z, i/ie function z^ a (z) has a 
unique pole in the strip S n . This pole is simple, and lies in the interior of S n . 

For the proof we need the following lemma. We will use the definition (|4.4|) of 
r(s, z) for general z € C \ M<o, the integral being over the infinite ray u £ z + M>o- 

Lemma 5. For any s € [1,2] and any z = —x + iy G C, satisfying either ^tt < \y\ < 
\x, jtt < \y\ < x or [x > and \y\ > ir], we have 

(A.5) |r(-«,*)| < s'^pV. 

Proof Take s and z = —x + zy satisfying the assumptions. By symmetry, we may 
assume y > 0. We may deform the contour of integration in (|4.4p to be the ray 
{z + f(l + fcz) : t > 0}, where k is any fixed non-negative number. This ray intersects 
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the imaginary axis at (y + kx)i, and thus \u\ > (y + kx)(l + k 2 )" 1 ! 2 holds for every 
point u on the ray, and 

Applying this with k = 1, we see that (|A.5|) holds whenever < ^r^- But 

> 1 for all non-zero z and thus the inequality holds for all s G [1,2] if and only 
if it holds for s = 2, i.e. if and only if r^jp < g- However, it is easily verified that 

CE I fit ^ Q 

^ x+ yyi is a decreasing function of x > for any fixed y > 0. Hence, if x > y > |7r, 
then f^±4 < i < i < i. similarly, if x > 2y > tt, then < —V < ± 

and if y > tt and x > |y, then ^_^ 3 < gjf^ < jtj|^ < g. To treat the remaining 
case, when y > tt and < x < |y, we apply (|A.6p with k = 0; from this we see 
that ()A.5p holds whenever s(|z|/t/) s < y. However, if y > tt and < x < |y, then 
s{\z\/y) s < 2{\z\/y) 2 < ^ < tt < y, and we are done. □ 

We also record the following bound, which follows from (|A.6|) with k = 1: 

Lemma 6. T/ie bound \T(—s,z)\ <C |,z|~ s_1 e _3f ^ ZioZrfs uniformly for all s G [1,2] 
and all z £ C with Ikz < 0, ^ 0. 

Proof of Proposition^ Let r/a(z) = z»<I> a (2) = e lz — iz Jq 1 e lzt t~a dt and note that 
77 a is an entire function. By (|A.1[) . our task is to prove that for each n, rf a (z) has a 
unique zero in S n , which is simple and lies in the interior of S n . Using (|4.21|) and 
applying the recursion formula T(s, z) = e~ z z s ~ 1 + (s — l)T(s — 1, z) twice, we find 
that for z with > 0, we have 

w\ = (— iz)°T(l — a -1 ) 
(A. 7) r] a (z) = w\ + u>2 + W3 with < W2 = a~ l (— iz)~ l e %z 

w 3 = -2d£(-iz) l «T(-l - a-\ -iz), 

wherein {—iz)a = exp(^ log(— iz)) with the principal branch of the logarithm; — tt < 
9flog(-w)<0. 

Let n G Z + and z = x — iy G r ra . We wish to apply Lemma [5] with s = 1 + a -1 
and with — iz in place of z. In order to justify this application, we have to check 
that either x > 7r, y > x > §7r or y > 2x > tt; this is clear if n > 2, since 
then x > 7r, and if n = 1, then the claim follows using (|AT4l) . tan(ivr) = 1 and 
tan(g7r) > 2. The conclusion from Lemma [5] is that \w%\ < \w^\ holds in (|A.7|) . We 
also note that &rg(wi/w2) G (1 + a _1 )(— ^tt + axg(z)) — x + 2tt%, and by (|A.4j) . we 
have x G ((2n — |)7r, (2n — \)tt\ and arg(z) = — (n — \)tt + \x G ( — -^7r, 0]; together 
these imply that arg(— ize~ lz ui\) lies in [— io _1 7r, (| — a _1 )7r) C (— |vr, ^7r), i.e. that 
lk{w\/w2) > 0. Moreover, \w3\ < \w2\ forces $t((u)2+W3) /W2) > 0; hence we conclude 
that $i((wi + W2 + Ws)/w2) > 0, i.e. that 

(A.8) K(-^e -fe r? a (^)) > for all z G F n . 

This shows that r] a {z) has no zeros along T n , and also gives a precise control on the 
variation of argr/ a (z) along T n . 
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Next, from ()A.7|) and Lemma [6l we see that for z = x — iy with y large and x > 
bounded, we have r] a (z) = w\ + W2 + u>3 = u^(l + 0(y~ 1 )), and thus argr/ a (z) £ 
7r + x + 0(y _1 ) + 27rZ. Also note that 3fy a (z) > for all z > 0, since K$ a (z) > for 
all z > (as noted previously) and r] a (0) = 1. Using these facts together with (|A.8|) 
(applied both for n and n + 1), we conclude that for any n G Z + and any sufficiently 
large Y > (depending on both a and n), argr/ a (z) increases by 2tt as z travels 
around the boundary of S n n {9z > — Y} in the positive direction. Hence, by the 
argument principle, r] a (z) has a unique simple zero in the interior of S n . Using the 
symmetry rj a (—z) = r) a (z), one proves the same fact also for So and any S n , n < 0. 
This completes the proof of the proposition. □ 

From now on, we write Cn = Cn( a ) f° r the unique pole of z^ a (z) in S n (n £ Z). 
By symmetry we have Q- n = — Q n for all n, and in particular £0 lies on the negative 
imaginary axis. Figure [2] below shows the curves traced by Qq, . . . , £4 as a varies. 

The next lemma gives an asymptotic formula for Q n (n > 0) with an error which 
is small whenever at least one of n, a and (a — l)^ 1 is large. 

Lemma 7. We have, uniformly over all a > 1 and all n £ Z + , 

(A.9) (n = (2n - a- 1 )^ + (1 + a" 1 ) arctan(^) - iY n + 0[ - jrr ] , 

V Y. n J \n + log|r(-a- i )| / 

where Y n equals the unique root y > of the equation 

(A.10) y-i(l + a- 1 )log((2vrn) 2 +y 2 ) =log|r(-a- 1 )|. 

(Regarding the error term in (|A.9p . we remark that ct 1 ) | > 3, and thus that 
log|r(-a _1 )| > 1, for all a > 1.) 

Proof. Using (|A.7j) and LemmaEJ together with the fact that ?? a (Cn) = 0, we get 

(A.ll) r(-a- 1 ) = (-< n )- 1 -^e^(l + 0(|Cn|- 1 )), Va>l, neZ+. 

Writing Q n = x n — iy n (x n ,y n > 0) and taking absolute values in (jA.lip . we get 
(A.12) \r(-a- r )\ = [Cr^-e^Cl + OdCnl" 1 )). 

Now, using the facts that \C, n \ > x n > (2n — |)7r S> n and |r(— a _1 )| — > 00 as a — > 1 + 
or a — > 00, we conclude that y n must be large whenever at least one oi n, a and 
(a — is large; and due to the form of the error term in ()A.9p . we may without 
loss of generality restrict to the case when this holds. Note that also |£ n | must be 
large, since \( n \ > y n . 

In more precise terms, we have, considering the logarithm of equation ()A.12p . 

(A.13) y n = |(1 + a" 1 ) log(x 2 n + y 2 n ) + ^^(-a" 1 )] + OdCnl" 1 ). 

In particular, using (2n — |)7r < x n < (2n + |)7r and also log(x 2 l + y 2 l ) < ^y n + 21ogn 
(which holds since y n is large), we conclude that 

(A. 14) y n x logn + log a - 1 )|; and thus \C, n \ X x n + y n X n + log|r(-a _1 )|. 

(Note: "X" means "both <C and >".) Now x 2 n + yl = ((2vrn) 2 + y 2 n )(l + 0{\C, n \~ 1 )), 
and thus, in (|A.13|) . we may replace "log(x 2 + y 2 )" by "log((27rn) 2 + y 2 )"; the error 
from this operation is subsumed in the error term 0(|C„| _1 ). We also note that the 
expression in the left-hand side of (|A.10p is an increasing function of y > 0, which 
is negative for small y and the derivative of which lies in the interval (1 — (27r) _1 , 1], 
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for all y > 0. It follows that Y n (in the statement of the lemma) is well-defined, and 
also that 

(A.15) y n = Y n + 0(\( n \- 1 ). 

Next, taking the argument of both sides of (|A.11[) . we get 

(A.16) x n = (1 - a" 1 )! + (1 + a' 1 ) arg(C„) + 2/cvr + OGCnT 1 ) for some k G Z, 

where -f < arg(C„) < 0. Clearly {2k-l)n -O^- 1 ) < x n < {2k+ ^Tr+OdCnl -1 ), 
and in fact, since — arg(£ n ) ^> y n |Cn| _1 and y n is large, we even have x n < (2k + |)7r. 
But also x n > (2n — |)7r; hence k >n. On the other hand, since £ n lies to the left of 
the curve r n+ i, we have x n < 2arg(£„) + (2n + |)7r, and using this fact in ()A. 16|) . we 
get (1 — a -1 ) arg(Cn) > (2(fc — n) — 1 — ^a~ 1 )'ir — 0(|Cn| -1 )- This forces k < n, since 
arg(£ n ) < and is large. Hence we have proved that k = n. Finally, using ()A.15p 
and (2n-|)vr < x n < (2n+|)vr, we get |arg( Cn)+a rctan(^)| < Y n \C n \~ 2 <C |Cn| _1 - 
Now (pOll follows from (pOEJh (lATHl) and (00311 . □ 

We may also remark that F n , as defined in Lemma satisfies 



n + G 



(A.17) r. - G + I±2L (l + $±fc£) + <?' 2 ) + o 

with G = log |T(— a _1 )|. This is proved by direct substitution in (jA.lOp . using the 
properties of the left-hand side in (|A.10j) noted in the proof of Lemma [71 

We will now change the contour in (|A.3p , Let a > 1 be given, and fix r > 
sufficiently small so that ()A.3|) holds. For n G Z + and Y > 0, we let z n y be the 
unique point where T n intersects {Sz = —Y}, and let C n y be the contour going 
from — oo to — (2n — ^)-7r along R, then along T_ n := {— z : z G T n } to — z nj y, 
then along {Qz = —Y} to z nj y, further along T n to [2n — h)ir, and finally along R 
to +oo. By the residue theorem and Proposition [3j for every n G Z + there is some 
Yo = Yo(a, n) > such that for F > Yq, we have 

(A.18) / V a (z)dz = ^- [ * a (z)dz- ]T Res z=Cm * a (z). 



m=l— n 



Now let wi,W2,W3 be as in (|A.7p , By Lemma [6] there is an N = N(a) G Z + such 
that |u>3| < § 1^2 1 f° r an -2 ^ r n , n > N. Using also !K(u>i/u>2) > for all z G T n , 
we get |?7a(z)| = \w± + W2 + > ^-\w\\ and thus \^> a (z)\ <C ti" 1_ «e S2 for all 
n > N and z G T n . Also, for any fixed a and n, we have |^a(^)| ^> F _1 e y for all 
z G C n y n {9=z = —Y} (cf. Lemma E] and (|A,7|> ): thus \^ a (z)\ < e" 2y for these z. 
The above bounds imply limn^oo (limy^oo f c |^ a (z)| |dz|) = 0, and so 



(A.19) 



P(<7{t, } > c) = — / * (z)cte = - hm £ Res 2=Cm * a (z) = £ 

*— ' T nrn — — v> <r> — — rv 



ae" 2 ^. 



Here the last equality follows from an easy calculation using (|4.9p and (jA.ip , noticing 
that the sum is absolutely convergent, since, by Lemma [JJ and (|A.17p . we have 

(A.20) |ae- 2lC "| < ae- 2G (|n| + G)~ 2(1+ ^, Va > 1, n E Z \ {0}. 
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One also checks that the formula ()A.19p may be differentiated termwise with respect 
to a, yielding 

oo , 

(A.21) /(c) = 2 e ~ 2Kn { 2ai {Ta Cn )~ l 

n=— oo 

for the density function (cf. (|5.ip ). 

For c not too large, the formula (|A.21|) can be used to compute /(c) numerically 
to a decent precision. We have implemented this in |25[ numdensity.mpl]. Our 
experiments indicate that for any given a > 1 (a = 2c) and n 6 Z + , the asymptotic 
formula in Lemma is sufficiently accurate so that it can be used as the initial value 
in the Newton iteration algorithm solving for $ a (z) = 0, with rapid convergence. 
Also, is computed using 

£c-^.-~(£..w)^ 

= -a^d^e-^+^fr'i-a- 1 ) - f°° e^u' 1 '^ (log u) du\ 

which most often can be evaluated very quickly via repeated integration by parts; 
in the remaining cases we use numerical integration. 

The data for the graph in Figure [T] can be found in [25, density.dat]; it was 
assembled by computing /(c) (c = ^a) for a = 1 + j^k, k = 1, 2, ... , 400. For each 
a- value we truncated the sum in (|A.21j) at |n| < 400 (using also the obvious n O —n 
symmetry). It turns out that the terms in (|A.21j) decay roughly as n~ 2 ^ l+ ^ as 
n — > oo (cf. ([A~20"1) h In particular we have slower convergence for larger a and this is 
seen in the computations: Our numerics indicate that we obtain the first few /(c)- 
values to within an absolute error < 10 , whereas for a near 5 (where /(c) ~ 0.05) 
the error is < 10 -6 . Of course the precision can be improved by including more 
terms in (|A.2ip . again cf. [251 numdensity.mpl]. 
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